Straining graphene by chemical vapor deposition growth on copper 
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Strain can be used as an alternate way to tune the electronic properties of graphene. Here we 
demonstrate that it is possible to tune the uniform strain of graphene simply by changing the 
chemical vapor deposition growth temperature of graphene on copper. Due to the cooling of the 
graphene on copper system, we can induce a uniform compressive strain on graphene. The strain is 
analyzed by Raman spectroscopy, where a shift in the 2D peak is observed and compared to our ab 
initio calculations of the graphene on copper system as a function of strain. 
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The discoveryof the electric field effect in a few layers 
of graphene 'l', '21 led to a genuine explosion of work on 
graphene with the subsequent observation of the anoma- 
lous quantum Hall effect in graphene monolayers 0, 0] 
following its theoretical prediction [s^. The most impor- 
tant property for electronic device applications is the pos- 
sibility to tune the resistance and electronic density of the 
graphene sheet. This is achieved by using the substrate 
directly as a backgate [l|, H, |3| or by processing a top gate 
i- 

More recently, it was suggested theoretically that it is 
possible to tune the electronic density of graphene and 
create confinement @, as well as to observe a zero-field 
quantum Hall effect by using strain engineering f^. In 
general, strain as a way to enhance electronic and opti- 
cal properties in semiconductors has been very successful 
Q . Mechanically, graphene turns out to be of very high 
strength |9i and is stretchable up to 30% in some cases 

While exfoliated graphene bears very high mobilities 
when suspended [111] , the size is limited to a few microns. 
Large scale epitaxial graphene on SiC is an interesting 
avenue but suffers from high substrate costs and difficulty 
to transfer Alternatively, large scale graphene can 

be synthesized by chemical vapor deposition (CVD) using 
a polycrystalline metal foil as a catalyst [13, ■ 

In this letter, we demonstrate that it is possible to 
strain graphene intrinsically and uniformly during the 
graphene growth process. This leads to interesting new 
properties, such as a change in Fermi energy and Fermi 
velocity in graphene. This opens new doors to the possi- 
bility of engineering devices, where the electronic prop- 
erties are tuned by strain. We consider the situation 
depicted in figure [U where a graphene monolayer is sit- 
ting on a Cu surface. For our ab initio calculations we 
consider a Cu(lll) surface, where the graphene and the 



FIG. 1: View of graphene on a Cu(lll) surface. The grey 
and dark orange spheres represent the C and Cu atoms, re- 
spectively. The blue parallelogram delimits the area of the 
supercell used for the first principles calculations. The letters 
A, B and C depict the different Cu(lll) sites. The minimal 
energy configuration of graphene on copper corresponds to 
one sub-lattice carbon sitting at the top-site (A-site) and the 
other sub-lattice carbon residing at the hollow-site (C-site). 



Cu share the same in-plane lattice constant We as- 
sume that the system has no structural defects, meaning 
that the graphene and the Cu match perfectly at the in- 
terface. In this case, the most stable configuration of 
graphene on Cu(lll) corresponds to one carbon sitting 
at the top-site (A-site) and the other carbon located at 



the hollow-site (C-site) (15|,|16|. The same configuration 
also applies to graphene on other metallic surfaces jT?}. 
To induce strain, we uniformly expand or compress the 
in-plane lattice constant of the whole graphene/Cu sys- 
tem. 

Zero strain corresponds to a C-C bond length of 
aeq=1.415A, which was found to minimize the total en- 
ergy of an isolated graphene sheet using density func- 
tional theory. For every a, we (i) perform a structural 
relaxation of the system to ensure that the forces acting 
on each atom is less than 0.01 eV/A, (ii) obtain the self- 
consistent electronic density and Hamiltonian and (iii) 
calculate the band structure of the hybrid graphene/Cu 
system along the high-symmetry points F, M and K. 

The most striking features induced by the strain are 
the change in Fermi velocity close to the Dirac point and 
more notably a shift in the Fermi energy with respect to 
the Dirac point as a function of strain as shown in figure 
[5]b). This is very interesting, since this allows for the 
tuning of the electronic density in graphene by strain. 

Experimentally, we expect to be able to strain 
graphene on Cu simply by heating the combined sys- 
tem, since the thermal expansion of Cu is different from 
that of graphene. While graphene has a small nega- 
tive thermal expansion coefficient a 18|, for Cu we have 
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FIG. 2; a) Band structure of graphene on Cu close to the K 
point for two different values of strain, i.e., 0% (solid line) and 
-6.7% (dashed line), b) depicts the change of Ef (relative to 
the Dirac point), 5k, 5q and KK' as a function of strain, c) 
Partial view of the Brillouin zone. The arrow labeled Sk^q 
indicates the phonon process at the origin of the unstrained 
2D Raman peak. The arrow labeled Sq^q corresponds to the 
wavenumber of the Raman excited electron of the unstrained 
lattice and 5q for the strained (compressed) lattice. The cir- 
cles show the constant energy surface of the Raman excited 
electron in the unstrained (solid line) and strained (dashed 
lined) lattice. 



a ~ 25.8 X 10-6 /°C [19.] at 1000°C. This would lead 
to a relative graphene compression of 0.5% from 1100 to 
900°C. 

We synthesized graphene monolayers by CVD of hy- 
drocarbons on 25 /Ltm-thick commercial polycrystalline 
copper foils. The Cu foil is first acid-treated for 10 mins 
using acetic acid and then washed thoroughly with de- 
ionized water. Graphene growth is realized in similar 
conditions to the recently reported ones [2^ [2lj , but us- 
ing a vertical quartz tube. Graphene is grown with tem- 
peratures ranging from 900°C to 1100°C in 0.5 Torr, with 
a 4 seem H2 flow and a 40 seem CH4 flow for 30 minutes. 
During the cooling process the methane flow is stopped 
while the hydrogen flow is kept on. 

Figure [3] illustrates the temperature effect on the syn- 
thesis of graphene. These scanning electron micro- 
scope (SEM) images show that the density and size of 
growth domains increase with temperature. At the low- 
est growth temperatures we did not observe full coverage, 
even for longer growth times. 

In order to characterize the graphene layer further, 
we used Raman spectroscopy. Raman spectra are mea- 
sured with a lOOx objective at 514nm, having a 1800 
grooves/mm grating and spectral resolution of about 2 
cm"^. The results are shown in figure ID There are two 
dominant peaks, labeled G and 2D. The G peak around 
1576 cm"^ arises from a first order G band phonon pro- 
cess. The energy corresponds to the phonon energy at the 





FIG. 3: SEM images of graphene on Cu grown for 30 minutes 
with different growth temperatures: (a) 900°C, (b) 950°C (c) 
975° C, and (d) 1050° C, respectively. Each image is approxi- 
mately 30/imx30/xm. 



F point, which is degenerate. Several authors have shown 
that the degeneracy is lifted by the application of uniax- 
ial strain [22], which leads to a splitting of the G peak. 
For uniform biaxial strain, there are no experiments on 
compressive strain in graphene, but we expect the strain 
to simply renormalize the phonon energy {to'^) following 
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1 - 2e7) |23|, |2J]. 7 ~ 1.24 is the Griineisen 



IS 



parameter in carbon nanotubes 25|, e = (a — aeq)/a, 
the strain and the factor 2 stems from the biaxial strain 
in both directions. 
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FIG. 4: Raman spectra for each growth temperature normal- 
ized to the 2D peak. Inset: Close-up of the 2D peaks showing 
the shifts of the peaks. The 975° C spectrum shows a region 
of possible bilayer. 

The other dominant peak is labeled 2D (also called 
G') and stems from a two phonon process. The two 
phonons of wavenumber Sk and -Sk scatter the electron 
with wavenumber Sq into the other valley (K to K') 
and back as illustrated in figure [2] c). The phonon en- 
ergy depends on the dispersion close to the Dirac points 
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(K or K') along the M direction. Hence Sk will de- 
pend on the electronic band structure, which is shown 
in figure [2] a) for two different values of strain. Us- 
ing the calculated electronic band structure we can ex- 
tract Sq in the M direction for our Raman laser energy 
Eiaser = 2.412 eV (514 nm). From a geometrical consid- 
eration, we have 5k = ^(KK')^ -|- 2 • KK' -Sq + A-Sq^ 
(see figure[2]c)), where KK'=K— K' is the distance in re- 
ciprocal space between the two valleys K and K'. KK' 
will simply be renormalized by a uniform strain leading 
to KK' = (KK')e,/(l + e), which is shown in figured 
b). For small values of e, the renormalization constant is 
approximately linear in strain. 5fc, on the other hand, is 
non-linear as it depends on 5q. 

In this region, the phonon band velocity can be pa- 
rameterized by ViTo — 5.47 x lO^'^vp 123[, where vp — 
lO^m/s is the Fermi velocity. In order to evaluate the 
change in phonon energy due to uniform strain for the 
2D peak we simply combine the change in wavenumber 
(5k from the linear dispersion with the renormalized to- 
tal phonon energy due to the Griineisen parameter as 
a;2^ cii [i^l^ + 2v,To{5k-5k^q)]{l-2e-i). The factor 2 in 
front of the band velocity term stems from the phonon 
pair involved in the process. For the equilibrium (un- 
strained) values of uJeq ' , we use those obtained for ex- 
foliated graphene at the same wavelength (514nm) [26| . 
i.e., = 2678cm~^ and = 1576cm~^. We can now 
evaluate the expected frequencies for the 2D and G peaks 
as a function of a uniform strain, which is shown in figure 
El The corresponding growth temperature dependence of 
the 2D peak data is then fitted to the 2D peak position 
obtained from the strain calculation, assuming a linear 
dependence. We also show the result (in dotted lines) for 
the 2D peak computed without the band structure ef- 
fect, which is then simply a linear function of strain. By 
including the band structure effect, the 2D peak depen- 
dence becomes non-linear with strain. This is in contrast 
to the G peak, which does not depend on the electronic 
band structure. 

The data fits nicely with the expected Raman shift 
when the difference in strain between growth tempera- 
tures of 900 and 1100°C is -0.5%, which is the expected 
value from the thermal expansion of copper for this tem- 
perature range. This clearly demonstrates that it is possi- 
ble to strain graphene uniformly, simply by changing the 
CVD growth temperature, which is one of the important 
results of this letter. This also shows that when graphene 
is grown at high temperatures it is under small compres- 
sive strain at room temperature (the Raman spectra are 
taken at room temperature). 

The peak labeled D in figure |4] is located at about half 
the frequency of the 2D peak, since it involves only a one 
phonon process. This process, however, is only activated 
in the presence of disorder |23l . |27[ . Raman spectroscopy 
can be employed to determine the number of graphene 
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FIG. 5: The Raman shift of the 2D and G peaks as a function 
of strain obtained from the ah initio results. The solid line and 
the dotted line for the 2D peak are obtained with and without 
considering the effect of the electronic band structure. The 
square symbols are the experimental data points obtained for 
different growth temperatures. The inset shows a zoom-in 
of the 2D peak data points and the line represents the best 
linear fit. 



layers [2g, which confirms that we have mainly single 
layered graphene as expected in CVD on Cu [20|. 

It is instructive to compare our results to those ob- 
tained for exfoliated graphene, where several groups have 
been able to detect uniaxial stain in exfoliated graphene 
flakes using Raman spectroscopy [2^. The reported 
shifts in the 2D peak vary from -7.8 cm~^/l% to -66 
cm~-^/l%. It was recently suggested that the large range 
of values could be due to the dependence on the polar- 
ization of the laser beam and orientation of the graphene 
lattice Our results, between and -1% strain, give 
values of -66 cm~^/l% (without band structure effect) 
and -79 cm~^/l% (with band structure effect), which 
is the higher bound of the reported values for uniaxial 
strain. This is consistent with expectations, since for bi- 
axial strain we expect approximately a doubling of the 
strain induced Raman shift. In a very different geometry, 
where stretched biaxial strain has been measured by de- 
positing a graphene flake on a shallow depression, a very 
high 2D peak shift of about 200 cm~^/l% was reported 
3ll |. In epitaxial graphene, shifts of the principal Ra- 
man peaks were observed to be dependent on the growth 
time for a fixed annealing temperature [32| . For CVD on 
nickel, no systematic shift with growth temperature was 
observed [33 1. 

The formed graphene/Cu heterostructure constitutes 
a very interesting system by itself. Not only can the 
amount of relative strain be tuned as shown above, lead- 
ing to changes in the Fermi energy and Fermi velocity. 
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but this system also opens the door to numerous appli- 
cations. For instance, the added graphene layer can en- 
hance electronic properties in Cu interconnects, for ex- 
ample, or greatly enhance thermal conductivities of thin 
Cu films. The tunability of the strain of the graphene 
layer, greatly increases this potential. The graphene/Cu 
heterostructure can also be selectively etched for partial 
or full transfer onto other substrates, which can lead to 
interesting transport properties such as a very sharp weak 
localization peak 34| . 
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